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We introduce and analyze several low-dimensional scattering systems that exhibit geometric phase 
phenomena. The systems are fully solvable and we compare exact solutions of them with those ob- 
tained in a Born-Oppenheimer projection approximation. We illustrate how geometric magnetism 
manifests in them, and explore the relationship between solutions obtained in the diabatic and adi- 
abatic pictures. We provide an example, involving a neutral atom dressed by an external field, in 
which the system mimics the behavior of a charged particle that interacts with, and is scattered by, 
a ferromagnetic material. We also introduce a similar system that exhibits Aharonov-Bohm scatter- 
ing. We propose some practical applications. We provide a theoretical approach that underscores 
universality in the appearance of geometric gauge forces. We do not insist on degeneracies in the 
adiabatic Hamiltonian, and we posit that the emergence of geometric gauge forces is a consequence 
of symmetry breaking in the latter. 
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I. INTRODUCTION 

The Born-Oppenheimer (BO) approximation allows the replacement of a complex many-body problem by a mean 
field description in which the motion of atoms can be accurately described. The dynamics of the latter are governed 
by an effective force F generated by a scalar potential that represents the collective motion of the fast degrees of 
freedom, i.e. the electrons. According to the Hellmann-Feynman theorem^ 

F^-Ve{R) (1) 

where e is the Born-Oppenheimer energy of the electronic state and R typically describes the geometry of the atomic 
ensemble. Molecular spectroscopy, atomic collision phenomena, as well as molecular dynamics simulations are just 
some of the applications that derive their utility from the implications of Eq. Q. In its modern expression [5H3] BO 
theory also allows for an additional contribution to Eq. ([T]), that of an induced effective "magnetic" force [1H7] 

„ V „ ^ V 

F^ = -xB-Bx- 
2 2 

B = V xhA (2) 

where v is the atom velocity operator [4| and A is a vector potential. This effective force is sometimes called geometric 
magnetism, a term coined by M.V. Berry, and is the nomenclature we adopt here. 

Advances [HHI3] in the manipulation of cold atomic matter has allowed researchers to "dress" atoms via the appli- 
cation of laser fields in order to induce gauge potentials that manifest as effective Lorentz forces on the atoms. This 
phenomenon has been called, "artificial" ^2j or "synthetic" [lOi magnetism and in a laboratory realization of it [8,, a 
BEG cloud was observed to undergo cyclotron motion in the same manner charged particles behave in a magnetic 
field. The effect has important implications for the control of atomic matter, and offers a novel route to quantum 
magnetism [m [M] . 

In this paper we address the question; what is the relationship between geometric and induced synthetic magnetism? 
The appearance of vector gauge potentials in molecular physics was first demonstrated by Mead and Thrular[T5j in 
their study of pathologies associated with conical intersections (degeneracies of distinct BO molecular energies at 
an isolated point Rc)- At a conical intersection, the phase of the BO wavefunction is undefined and, in order to 
avoid multi-valuedness, a compensating phase factor must accompany the nuclear wavefunction. That procedure 
is equivalent to the requirement that a non-trivial vector potential is minimally coupled to the nuclear motion JS]. 
The Mead-Thrular gauge potentials do not exhibit geometric magnetism but since Berry's vector potential[5! I16j has 
non-vanishing curl, systems exhibiting the latter do. In the Mead-Thrular treatment degeneracy plays a central role 
as it does in Berry's IIamiltonian[T2]. Degeneracy is also an essential ingredient in the non-Abelian generalization of 
Berry's phase[T2. Thus the appearance of geometric phases has historically been associated with systems in which 
degeneracies manifest. This does not appear to be the case for artificial magnetism[T2] and so it is worthwhile to 
further investigate the relationship between the two phenomena. 

In applications [HHSnj it is often necessary to go beyond BO theory and exploit a multi-channel generalization of 
it, the Born-Huang expansion or the method of perturbed stationary states (PSS)[2T]. In atomic collision theory the 
latter lead to equations that have the formjl] 

(V -iAfF + V F ^ EF (3) 

2m 

where -F is a multi-channel amplitude expressed as a column matrix, A is a matrix-valued vector potential, and V 
is a diagonal scalar-valued matrix whose elements consist of the BO eigenvalues associated with each channel and E 
is the total energy eigenvalue. If a complete set of channel states are included in the PSS expansion then Eq. (|3| is 
exact. We distinguish two cases in applications of Eq. ([3|. In the first, matrix V is degenerate at all R and a gauge 
transformation can be found[4l|6] so that, in the new gauge, A = 0. The representation in which A = is sometimes 
called the diabatic picture[ll|^, whereas Eq. ^ describes dynamics in the adiabatic picture. The second case, which 
we deal with exclusively in this paper, V is non-degenerate and Eq. ^ exhibits gauge covariance if the vector valued 
potential A is taken as the spatial component of a 3-j-l gauge field that also contains a temporal component|3]. In 
applications [T8'-'20^ one typically truncates the PSS equations to a finite number of channels (most realistic systems 
involve an infinite number of channels) and in that case one cannot, in general, find a transformation into a gauge 
where A = 0. As a consequence effective gauge forces [H [TS] manifest. Does this imply that the existence of the latter 
are the artifacts of an approximation procedure? What is the relationship between geometric forces and the fictitious 
forces that arise in a classical description where non-inertial frames are employed? In order to address these, and 
other questions, we introduce and analyze several low dimensional but solvable models that exhibit geometric phase 
phenomena in the BO approximation. In the systems discussed here the "fast" degrees of freedom span a finite Hilbert 
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space of dimensions d = 2, 4. Similar models have been discussed previously [23'-'5S] but here we focus on systems 
that do not exhibit energy degeneracies in the adiabatic approximation. We concentrate our efforts on fully quantal 
treatments of Eq. (|3| and we limit our discussion to scattering systems and their solution. 

In section I we provide a general theoretical framework for the subsequent discussions. In section II we introduce 
a ID scattering system in which geometric phase phenomena arise and which fully quantal analytic solutions are 
available. In section III we introduce a 2D scattering system that exhibits the phenomenon of synthetic, or geometric, 
magnetism. We compare the scattering solutions of it with that obtained for a charged particle scattered by a 
ferromagnetic slab. In the remaining sections wc show how these systems are related to ones in which geometric forces 
arise due to inter-particle interactions. An example involving two, interacting, spin-1/2 particles is presented. 

Vector valued quantities will be shown in boldface, we will announce matrix-valued quantities in the text and use 
ordinary Roman typeface to represent them. However, in certain circumstance, where there might be a possibility for 
confusion, we will underline a Roman letter to stress its matrix nature. 



II. THEORY 



Consider the Hamiltonian 

H = -^V^ + Had{R) (4) 

where Had{R) is the Hamiltonian spanning an rt-dimensional Hilbert space /i„ that represents the internal degrees 
of freedom of a quantum system of mass m. It is parameterized by the eigenvalues of the quantum variable R. In 
molecular physics Had is called the adiabatic Hamiltonian and represents the total kinetic, electrostatic, and magnetic 
interactions among the electrons. In that case the internal space is spanned by the Born-Oppenheimer (BO) 
eigenstates of Had- Additional realizations described by Eq. Q could be an atom, molecule or spin-n system in an 
external field that is modulated by the values of R. If we assume that the basis for /i„ is discrete, we can express 

Had{R) = U{R)HBoU\R) (5) 

where Hbo is a diagonal matrix whose entries we shall label ei, 62, ...e„. In general they are parameterized by R (more 
precisely, by the eigenvalues of the quantum variable R). In molecular physics they are called the BO eigenvalues 
for the electronic Hamiltonian Had- U{R) is a unitary operator acting on the internal states and, in general, is also 
parameterized by R. Expressed as a matrix representation in the BO basis, [/ is an n x n unitary matrix. 

Any well behaved Hermitian operator Had can be written in the form given by Eq. ([5|, however we require the 
additional condition that U must be single-valued in the parameter space spanned by R. We now posit the following. 
Any such unitary matrix can be written in the form, 



U'' {R) = V exp(i j dR' - A{R!)) 



(6) 



where A is a vector valued n x n matrix (i.e. a non-Abelian gauge potential), and V represents a path-ordering 
operator. Path-ordered integrals along curve C in parameter space R are defined as follows: Consider a path, C, 
defined by the set Ri,R2, - Rn where Ri = R{ti),R2 = R{t2)---Rn = R{tn), Rn = R, and ii < t2 < defines a 
trajectory that maps out the path. We define 



z-*" dR{t) 



Vexpii [ dR ■ A{R')) = Texp{i [ dt ■ 
Jc Jti 



dt 



A{t)) (7) 



where we assumed that the path is sufficiently smooth so that ^^^j^ is well defined on it and T is the Dyson time- 
ordering operator. If C traces a closed path Eq. ^ represents the Wilson loop integrali^fij which we require to vanish. 
To that end, we demand that [4] [27] 

F^,, = duA, - d,A^ - i[A^, A,] = 0. (8) 

We make the assumption, if A satisfies condition Eq. ([s]) and is not singular, then the value for path integral Eq. Q 
is independent of C. With this working assumption, we are allowed to take the gradient of Eq. Q 

{yU^)U^iA (9) 
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i'VU^)U = -U^VU = lA (10) 

and 

A = iU^VU. (11) 

This proves our assertion, given the unitary operator U{R) we can always find a gauge potential given by expression 
(11), that satisfies identity ([6|. An illustration is the direct product of an /i„ = 2 dimensional Hilbert space of a 



spin-1/2 system and a ID Hilbert space for quantum variable i2, e.g. R <Z x and —oo < x < oo. We define the gauge 
field 

1 -Aq 



where Aq,Ai are real constants "55". It is evident that condition Eq. Q is satisfied by this gauge potential and we 
can apply Eq. ([6| to construct 

/ cos {Ax) - i-Ao^Ax _rMj^Ax 

\ ^lAx^fAx COS (Aa;) + 
A = /V7V. (13) 
We note that U is single- valued and U{x = Q) is the identity operator. Suppose the internal Hamiltonian is given by 

Had{x) - UqHboUI (14) 

where Hgo is an arbitrary diagonal matrix. Previous studies e.g. [151 123j have focused on systems in which Hgo is 
degenerate for some value of the parameter i2c, i.e. ei(i?c) — 62 (-Rc)- Here we do not allow any crossings, indeed 
we take the BO eigenergies to be constant throughout the domain of x. We let 

^""(oi) ^"^^ 

Uq^U X > (15) 



where U is given by Eq. (13 1. Choosing 



i^BO = ( j , (16) 
where A > is a constant, we obtain the set of coupled Schrodinger equations for the spinor eigenstates of Hamiltonian 

Z"-|j^Z+|ji^Z = o (17) 



where 



and 



V^UHboU^ x>0. (19) 
We seek scattering solutions to Eq. \\7\ and if we set E = —A + where < A, the excited internal BO 



channel is closed. Instead of solving Eqs. (17) directly, it is convenient to transform to the adiabatic picture whose 
wave function is. 

Fad = UoF (20) 
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and satisfies 15 



(V - lAfF,^ - ^ HboF,^ 



2m 
7i2 



(21) 



where A is the gauge potential 



A: 
A : 



^0 Ai 

Ai -Ao 

a; < 0. 



a; > 



(22) 



In the adiabatic picture we have replaced the ofF-diagonal potential matrix V with the diagonal BO matrix, for 
the price of gauge potentials. At collision energies where some of the excited BO channels are closed, a common 
approximation in molecular and collision physics is to project the system of coupled equations (which are exact) 
onto the open sector. This is called the Born-Oppenheimer (BO) or, if several channels are open, the perturbed 
stationary states (PSS) approximation [U [31]. In our example this approximation leads to the following single channel 
equation [4], 



{d + iAofFo + k'^Fo-hFo = Q x>Q 
d^Fo + k^Fo = a; < 



(23) 



The induced potential b is closely related to the so-called "B" term, or adiabatic, correction [29H3T]. A gauge transfor- 
mation Fo — >■ exp(zAoa;)Fo allows us to "eliminate" the derivative coupling in Eq. (23 1 and we can proceed to solve 
the scattering problem. Imposing the boundary condition that the open-channel wavefunction vanishes at x = L >> 
(which is equivalent to the placement of an impenetrable barrier at L), we obtain, in the asymptotic region for x < 0, 



Fq{x) = exp{ik{x — L)) + Rexp{—ik{x — L)) 

2k 

R = exp(-2ifcL) I -1 + 



R = exp(-2iA:L) -1 



k + i^k^-Aj cot (l - A?) 
2k 

k + i^A\ - fc2 coth {l^A\ ~ fc2 



fc> lAil 



fc< lAil 



(24) 



We can also solve the fully coupled Eqs. (21) analytically, as outlined in the appendix, without resorting to the 



Born-Oppenheimer approximation. That solution, in the low energy limit, leads to 



i? = -1 + 2ifc L - 



tanh(AiL) 



0(fc2) 



(25) 



where we have kept only the lowest order terms, as fc — > 0, in an effective range expansion for R. This result is in 
harmony with that obtained in an effective range expansion of expression (24) obtained using the BO approximation. 
In summary, 

• At collision energies where the excited internal state is closed and the energy defect between the open and closed 
channels is large, the lowest order term in an effective range expansion (i.e. in the cold collision energy regime) 
of the Born-Oppenheimer expression for the reflection coefhcient is identical to that obtained by solution of the 
fully coupled equations. It depends on the value of the off-diagonal component of the induced gauge potential 
( 12 ) despite the fact that the pure gauge condition. Eq. ([s]) is satisfied. 



• In the BO approximation the off-diagonal gauge coupling leads to a higher order induced scalar potential, given 
by the expression for h in Eq. ( Al ) . In this approximation the diagonal component Aq can simply be "gauged" 



away and physical quantities, such as the scattering length, are independent of it. 

• At collision regimes where the BO approximation is no longer appropriate [32], (e.g. when A ^ i?), the diagonal 
components do affect scattering properties. This may seem counter-intuitive since we stated that a gauge 
transformation can be employed to transform these potentials away. One way of proceeding beyond BO theory 
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is to describe the couphng between the open and closed channel through the introduction of a non-local optical 
potential's]. In such a treatment Eq. (23 1 is replaced by [34] 



2m 

{d + iA^fFo + k'Fo + VoptFo = x>0 
VoptFo= I dx'Voptix,x')Foix') 



(26) 



where Vopt is a nonlocal potential. Now the transformation Fo — > exp{iAo x)Fo does eliminate the derivative 
coupling in Eq. (A7| but, since Vopt is non-local, it results in 



dx' Vopt{x,x')Fo{x') I dx' Vopt{x,x')e-x:p{iAQx')Fo(x') 



and it cannot simply be "gauged" away as in the BO approximation. 

• If the collision energy is sufficiently large so that i? >> A and both channels are open the resulting scattering 
properties become independent of both Aq and Ai . In the case where the internal state energy defect A can be 
neglected, we can treat A as a pure gauge, i.e. it does not affect the scattering properties of the system. 



III. GEOMETRIC MAGNETISM 



The model introduced above has interesting and salient properties but it lacks the feature of geometric magnetism. 
Because the diagonal components can be gauged away in the adiabatic, or BO, limit effective Lorentz forces do 
not manifest. Can we construct a gauge potential (again here we limit ourselves to the /12 or spin-1/2 internal space) 
that satisfies the integrability condition Eq. ([s]), but includes diagonal components that are non-trivial? The answer 
is yes but we need to extend our parameter space, at the least, into two dimensions i.e R G x,y. Consider 



U = cxp(-ia-30(x, y)) exp{-ia2i^ix, y)) exp(i(T30(a;, y)) 



(27) 



here ai are the Pauli spin matrices and (f){x,y) and Q{x,y) are arbitrary but single-valued functions in the xy-plane. 
We have chosen this form because its structure is similar to that which manifest in molecular physics'18-'20' . By 
construction U is single valued and the gauge potentials derived from it satisfy integrability condition Eq. (|8| . Using 
prescription (111 we find 



A = 



-2sin2(r2)V0 
exp(2i0) {iVn - sin(2f^)V(/ 



exp(- 



ivn) 



'2i(j}) (sin(2f7)V(; 

There is freedom of choice in VL, 0, as long as single-valuedness is enforced, but here let's introduce 

0(x,y) -LB02//2 



(28) 



^{x,y) = 9h{x)9h{L - x)Arcsin{yJx/ L) + 9h{x - i)^ 



(29) 



(A similar but not equivalent system was independently proposed in Ref. !l4j) where 9h is the Heaviside step function 
and Bq,L are positive constants. Inserting Eq. (29 1 into Eq. (28 1 we obtain, For a; < L, A = 0, in the region 
< a; < L, 



A 



( 

V 



,2i4 



-X Boj 
-sin(2f7)^j - 



-2i<f> 



sin(217)^fa3' 



xB^j 



and for x > L 



-BqL 
SoL 



(30) 



(31) 



We keep Hbo described in the previous paragraphs. In seeking scattering solutions for this system it is again 
convenient to proceed in the adiabatic picture. The adiabatic amplitudes obey equations (21 1 except that the gauge 
potentials are now given by Eq. (28 1. As above, the ground internal BO state is open and the excited state is closed 
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and it is appropriate, for sufficiently low collision energies and a large energy gap A, to proceed with a BO projection 
of Eqs. (21 1 into a single channel description. Hence, for < a; < L, we obtain 

d^Fo + {dy ~ iA^fFo + k^Fo - h{x) Fa = Q 
Bqx 

b{x) = ■ Am = [sn?{2n)V<l) ■V(j)+Vn- VVl] = 

= <L-x)Bl^^-^^. (32) 

This Schrodinger equation is equivalent to that obtained for a system in which a charged particle interacts with an 
effective magnetic field 

B = V X hAo = k HBq <x < L, (33) 

in addition to an induced scalar potential b{x). The proposed spin 1/2 system does not have an electromagnetic charge 
and so this effective magnetic field is a manifestation of what is, sometimes, called artificial or synthetic magnetism[Sl- 
IT^ [55] . It leads to effective Lorentz type forces acting on the dynamical system. The emergence of such forces was 
first demonstrated in molecular physics [IHB]. 



A. Non-singular model 

Whereas the above model describes quantum magnetism the induced scalar potential suffer singularities at x = 
and at a; = L and so we wish to introduce a similar model that does not posses these singularities. We now posit that 

(j>{x,y) = LBay/2 

n{x, y)^^(l + tanh[/3a;]) . (34) 

We then obtain for the vector potential, 

. _ / -Boisin(|7r(tanh(a;/3) + 1)) -^e-^^°^y cos {^ntanh{x 13)) \ 

^ \-^e'^°^ycos{lTrtanh{xl3)) BqL sin (i7r(tanh(a;/3) + 1)) J ^ ^ 

f -iie-'^°^Wsec/i2(x/3)\ 

[lie'B»^yn(3sech^ixl3) J' ^^^^ 



and 



In the limit a; — > cx3, 



- ( bIl ) ' (37) 



The diagonal components of the "magnetic" induction |-B| = fi?^ is given by 



whereas A^ — )• 0, and A — > as x — t- — oo 

leti 

h 

= =F 

where the total magnetic flux density is 



B{x) = t'-^ttISBoL sech^{(3x) cos (^- tanh(/3x) j (38) 



d,^/ dx^-^ = \B,\L (39) 



In the BO approximation, for the open channel, we obtain 

dlFo + {dy - lAofFo + k'Fo - h{x) F^ = 

Aq{x) ^ $ sin (^TT{tanh{x 13) + 1) 



b{x) ^J2^^k■ Aki = s\T?{2n)v<p ■V(j) + vn-vn = 

= — (2<I>2(cos(7rtanh(;3a;)) + 1) + Tr^/?^ sech'^ipx)) . (40) 
16 
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We assume the incoming wave is incident in the normal direction, i.e. Fq = Fo{x) and the effective equation for the 
BO amplitude is, 

dlFo - Alix) Fo + k'^Fo - b{x) F, = 0, (41) 

or 

dlFo + k^Fo~v,,ff{x)Fo^O 

Veffix) = Alix) + ^ (2$2(cos(7rtanh(/3x)) + 1) + ir^P^ sech'^{/3x)) . (42) 

In Figure [T] we plot the effective potential for a given /3 and $ = BqL. We note that transmission occurs when the 
collision energy E is greater than the height of the potential barrier $ as x — cx). In that diagram we superimpose 
the effective potential discussed in the appendix for the case of a ferromagnetic slab characterized by the same value 
of In the classical description the value $ determines whether the particle has sufficient velocity so that its Larmor 
radius is greater than the slab thickness and can thus penetrate the slab. Since A — >■ 0, as x — > — oo we impose the 





1 
1 

I. ± 


k 




1 i 






/ ./ 1 * 





FIG. 1. (Color online) Plot of the effective potential given by Eq. (42 1 (solid blue line). The dashed blue line represents the 



effective scalar potential h{x) given by Eq. (381. The dotted line represents the effective potential for the ferromagnetic slab 



discussed in the Appendix. The shaded region represent the slab. We have chosen the parameters so that both cases have 
identical values for the total integrated flux density defined by Eq. (39 1. 



incoming boundary condition 



Fq{x) — >■ exp(ifca;) + r exp(— ifcx) 



(43) 



where r is the reflection coefficient. If _B > $ transmission into the asymptotic region a; — >■ oo is allowed and, in this 
region, the vector potential has the form given by Eq. ( 37 ) . The outgoing asymptotic flux is given by (see appendix 
A) 



jy = --AyF:F,^-h\^^ 
ni m 

where we made use, for fc > |<E>|, of the scattering boundary conditions 

Fo{x) texp{i\/k'^ - ^"^x). 
The angle of deflection, with respect to the normal is then given by 

jy _ * 



tan( 



Jx 



V/c2 _ $2 



(44) 



(45) 



(46) 



and agrees with the result obtained in the classical description of scattering of a charged particle particle by a 
ferromagnetic slab (see Appendix). 
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B. Multichannel description 



The discussion presented above makes explicit use of the BO approximation, we now repeat the analysis using the 
full, and exact [36 , multichannel description. We need to solve the coupled equations, 

~2to 



F 
V 



V-iA] F + VF-EF = 

A 

-A 



(47) 



where A is given by Eqs. (351 and (36 1. The diabatic gauge amplitude G is related to F via gauge transformation 
Eq. (27) and it obeys 



2m 

w = uvu'f 



WG-EG=0 

Acos(2f](a;)) e-**?^A sin(217(2;)) 
e**?^Asin(2f7(a;)) -A cos(2f7(x)) 



(48) 



In Figure [T] we plot both the diagonal and off -diagonal components of the matrix-valued potential W using the 
expression for rt{x) given by Eq. (34) for y = 0. It is somewhat more convenient to define a new amplitude so that 




FIG. 2. (Color online) Schematic diagram of scattering setup described by Eqs. (48). (a) An incident wave in the open channel 
whose potential energy is given by the solid blue line. The horizontal axis represent the x coordinate, and for x < the channel 
described by the red solid line is closed. For a; > the transmitted wave is propagated on the potential surface given by the red 
solid line. The brown line represent the total collision energy of the system and the green dashed line represents the off-diagonal 
coupling between the two potential surfaces, (b) The same system now illustrated in the adiabatic picture. The blue line is 
the BO energy for the open channel and the red line represents the BO energy for the closed channel. 

G=UG' where 

U = exp(-i^a3) (49) 



and the equation for G' is identical to Eq. ( 48 ) except that W is replaced by 

(50) 



Asm(2i2(a;)) — A cos(2S2(a;)) 



Thus 



- — (V-iA] G' + W'G- EG' = 
2m ^ / 

A^iWVU = j \ ^ (51) 
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or, letting G' = ^ : 



+ 'I + 4>'-^' + I? i»'(2Si)/2 + |f - I? Asin(2Si)/, = 0, (52) 

With the ansatz fi{xy) = cxp{^i'^y)fi{x) and f2{xy) = cxp(— z^j/)/2(a;) we obtain, 



" <i>Vi - |f Acos(21])A + I? ^^/i - |J Asin(21])/2 = 

+ |J A cosmh + |J - |J A sin(2f7) A ^ 0. (53) 



In the asymptotic limit x — — oo, 17 — 0, and Eqs. (53) reduce to 

— -<^f,-—^f, + —Eh^Q 
d^f2 2m ^ . 2m 



A/2 + —Eh = Q (54) 



whose solutions are. 



G' - exp(-i-y) r exp(Ka:) \ 

— 2 \ exp(ifca;) + rii exp(— i/cx) / 



|j(A-i?) + $2 



We have chosen this boundary condition so that 



fcEE W^(£; + A). (55) 



G = C/^G' = exp(-i — cr )G = r exp(-i$y) ri2 exp(Kx) 
— — 2 ^ \ exp(?fcx) + rii exp(— zfca;) 



(56) 



describes an incoming, normally incident, plane wave in the open channel. In the region x — > oo, f2 — ^, and Eq. 
( 52 1 reduces to 

0-|?A/-t-/-O (.7, 

whose solutions can be written in terms of scattering boundary conditions 



e' = =p,-4,,(,';-;*v.)) 



or, 



k' = v/fc2 - $2 ^> = J^^A- E) (58) 



The conserved flux is given by 
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Inserting expression (56) into Eq. (60) we find that, in the asymptotic region x < 0, 

.U2^k + 2^k\rn\')^-(l-\rn\' 



Jx 



Jy 



ih 
2m' 



2m 



(^2iL exp(— 2Ka;)|ri2| 



(61) 



Like-wise in the region a; > 



Jy 



ih 

2rn 
ih 

2m 



[-2ik'\h2\ 

(2i<P\h2\ 



m 

-^\tl2\'. 

m 



We compare these solutions with those obtained in the BO approximation. According to Eq. (27 1, 

G = UF 



(62) 



(63) 



Fc 



exp(— i$ y) (^fi{x) cos fl — f2{x) sin fi^ 
/i {x) sin VL + f2 {x) cos 



Carrying out these transformations, using Eq. (56) and Eq. (|59[), we find as x — >■ — oo 



F = 



Fc 



and as 



Fo 



exp(— i$?/) ri2 exp(Kx) 
exp(ifcx) + rii exp(— ifcx) 



—ill exp(— Kx) exp(— i^y) 
ti2 exp(i-\/fc^ — '^^ x) 



(64) 



(65) 



(66) 



Projecting out the open channel amplitude (i.e. the Born-Oppenhcimer amplitude) we obtain 

Fo — >■ exp(ifcx) + rii exp(— ifcx) x — > — oo 
Fo ti2 exp{i\/k^~-W^x) x ^ oo. 

In Fig. (3) we present the result for the transmission coefficient 



ITI 



_ |il2|Vfc2-$2 



fc 



(67) 



(68) 



and compare it to the expression obtained using the BO approximation where ti2 is replaced by t given in Eq. (451 



The wavenumber fc is expressed in terms of the flux density $ and transmission is allowed for fc > $. In the coupled 
channel calculations we set the closed channel parameter k = fc, and so for larger fc the BO approximation becomes 
valid in the limit as k — >■ oo. It is clear from the figure that the predictions of the BO, or adiabatic, approximation are 
in very good agreement with the results obtained using the fully coupled calculations. As /3 get large the transmission 
coefficient is suppressed since the barrier height of the repulsive potential b{x), given by Eq. (401, becomes larger. 

We conclude that the result of the BO approximation, with the presence of a vector potential that describes 
geometric magnetism, does indeed provide an accurate account for the dynamics and is in harmony with the predictions 
of the exact treatment. 



C. Open channels 

According to the above discussion the amplitudes in the diabatic gauge are given by 



G= ( exp(-^$y)/i(x) I ^gg^ 
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Wavenumber k 



FIG. 3. (Color online) Transmission coefficient as a function of incident wavenumber k expressed in units of $ taken here to 
have the value of unity. The solid circles represent data obtained using the BO description, whereas the open circles denotes 
the transmission coefficient obtained using the coupled channel theory and given by Eq. ( 67 1 



where fi{x),f2{x) obey the coupled equations Eq. (52 1 and $ ^ BqL. The expression for the conserved flux is then 
given by 



TO V 



m 



Jy-^--\M' 



We now evaluate the integral 



J" EEE I ds- j 



(70) 



(71) 



for contour C shown in Figure |4j Consider the segments labeled I, III in the diagram. The contribution to Eq. ( 71 ) 
from these segments are 



h 

TO 



d2/'i>(l/i(a)p-|/i(-«)p), 



(72) 



like-wise for segment II we obtain 



TO 



(73) 



Since the integrand in Eq. (731 is independent of y this contribution is canceled by that along segment IV. Therefore, 
as \a\ — >■ oo we find 



J (J m TO 



(74) 



open. As X — > — oo Eq. (561 is replaced by 



where we made use of Eq. (|59|). We now consider the case where E > A + i.e. both channels in Eq. (|52|) are 



G 



exp(— i$2/) ri2 exp(— ifc2x) 
exp(ifca;) + rn exp(— ifcx) 



(75) 



where 



ft2^A/|j(i?-A)-a>2 



(76) 
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and which describes an incoming wave in channel 1 that can also reflect into channel 2. Similarly, in the region a; — > oo 
we require 



G 



exp(— ti2 exp(ifc'x) 
til exp(ifc2a;) 



(77) 



where 



k' = v/fc2~$2 



(78) 



We evaluate ds ■ j as above. Using Eqs. ( 70 ) we find that the contribution from segments II, IV cancel and that 
contributions from segments I, III give 



J" = / ds ■ j = —^2w 



l/i(a)p-|/i(-a)p 



-$2w 



TO 



\ti2\'-\ri2\' 



(79) 




FIG. 4. (Color online) Contour plot of flux given by Eq. (60 1. Shaded green area represents region where the magnetic 
induction B{x) given by Eq. (47 1 is non-negligible. Dashed blue line represents the contour C. Segment I, III extend from —w 
to w, and segments II, III range from —a to a. The collision energy is such that the upper adiabatic channel, shown in the 
diagram of Figure [l] is closed. 



To summarize, we found that 



2w^\ti2\ E<A 



TO 



1^121 



2to 



S > A 



2^2 



2m 



(80) 



D. Coupled equations in the adiabatic picture 



Before we discuss the significance of Eq. (80), it is useful to re-derive this result in the adiabatic representation. In 

conserved current 

V^t p _ V p 2iF^A f] 



it the amplitudes obey Eq. ( 47 1 and lead to the conserved current 

ih 



3 = 



2m 



where F is given by Eq. ( 64 ) . Consider first 



j = —(vF^ F-F'VF 
2m \ 



(81) 



(82) 
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We show below that 



J"; 



h 

( 

m 



ds- j = ^$2u;(|tnp-ki2|'l > A 



2m 



2m 



T^O E <A 



The contribution to the left-hand side of Eq. ( 83 ) from segment I is given by 

dyjy = - / dy{Im{Fo*dyFo) + Im{F/dyF,)) = / $|/i(-a)|^ 



where we have used Eq. (64 1 and the fact — > as a — > ~oo. Like-wise, for segment III we obtain 

h 



The contribution from segment II is 



dyjy = — / dy^\f2{a)\'. 
m 



dx{Im{Fo*d,Fo) + Im{F,*d,F,)), 



(83) 



(84) 



(85) 



(86) 



but since it is independent of y it cancels with the contr ibu tion arising from segment IV. Inse rting the limits for 
|/i(— a)p and 1/2(0)^, as \a\ 00, into Eq. (84| and Eq. (85) we arrive at our assertion Eq. (83). We now evaluate 

h 



m 

-^1 ds^ 
m Jc 



ds ■ F^AF = 



F,*F,Aii + F,*FoAi2 + Fo*F,A2i + F,* F,A 



22 



(87) 



where are the matrix elements given in Eq. (28 1. The contribution from segment II, 



i— / dx Fc*Fo exp(-i$w) - Fo* Fc exp(i$w) 
m J -a ^ 

dxd,n(^f*h-f;fi), 



is independent of y and it is canceled by the contribution from segment IV. Consider the contributions from segments 
I and III. First, we focus on terms that involve the off-diagonal elements, 

/ dysmi2n)(F,*F,expi-i<i>y) + Fo*F,exp{i<^>y)y (89) 

J — W 



m 2 



For segments I, III, 57 0, ^ respectively and according to Eq.(89l these contributions vanish. Finally, we consider 
the diagonal terms, 

h 



m 



$ / dysin^{n)(F,*F,-Fo*F, 

J —W 



whose non-vanishing contribution comes from segment III and has the value given by 



■Fa 



m 



dy[\h{a)\'-\f2{a)\'). 



(90) 



(91) 



Using Eqs. (64) and (69) we obtain 



^^=-$2u;(|ii2p-Kiip) E>A 
m \ / 

J-A = — $2u;|ii2p £;<A-|-- — . 
m Zm 



2m 



(92) 



Adding this contribution to that given by Eq. ( 83 1 we obtain 



F = F + Fa = -<^2w{ |ii2p - |ri2p 



m 

^ , ,9 . 

— $2w|ii2r -E<A-1 

m 2m 



E> A 



2m 



(93) 
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and which agrees with expression ( 80 1 obtained in the diabatic picture 



Equation (93) constitutes one of the main results in this paper. It expresses the scattering solutions in terms of the 
integral 

(94) 



T = J ds-j 



which is invariant with respect to the gauge or representation. That is, T has the same value in both the diabatic 
or adiabatic pictures. If the collision energy is such that the excited adiabatic channel is closed, then T does not 
vanish and its value is given by the second line in Eq. (93). Furthermore, in the adiabatic representation that value 
is determined from the contribution of the current that is proportional to A. Because this result is obtained in the, 
exact, coupled-channel framework it does not require the invocation of the BO approximation. However, the BO 



approximation provides an elegant and accurate description and the appearance of vector gauge potential in Eq. ( 40 ) 



is simply a consequence of the fact that the BO approximation is faithful to the requirement that J- does not have a 
null value. Indeed, the arbitrary removal of the vector potential in the BO approximation compromises its efficacy. 
It is important to note that the vector potential Eq. (40) and the resulting effective B field are source-less and 



there are no-singularities associated with it. Consider now the case where the BO approximation is no longer valid 
at higher collision energies where the excited adiabatic state is also open. In the high energy limit where E >> A we 
can ignore A (in lowest order) and it is evident that in this limit 



|ti2| \ri2\ 



(95) 



and — > 0. Thus, at higher collision energies the scattering properties are not affected by the presence of the vector 
gauge potential, i.e. its influence on the system is what one would expect from that of a pure gauge. This does not 
preclude the fact that when working in the adiabatic picture one cannot ignore the gauge potentials when solving for 
the multi-channel amplitudes. An interesting case concerns the energy regime where both channels are open, thus 
voiding the BO approximation, but the limit behavior, given in Eq. (95), has not been reached. That discussion will 
be deferred to future investigations by the author. 



In summary we have shown the necessity of including the non-trivial vector potential Eq, 
mation. The potential predicts a Lorentz force that determines the scattering properties of the proposed system 



(|40|) in the BO approxi- 
We 

illustrated how this system is an analog to that of a charged particle that is scattered by a ferromagnetic slab. We 
have solved the, exact, fully coupled equations in both the diabatic and adiabatic pictures and showed the fidelity of 
the BO approximation at low collision energies. At higher collision energies we find, in an appropriate limit, that the 
presence of the induced gauge potentials do not affect scattering properties. Thus the emergence of non-trivial gauge 
structure in the low energy limit is suggestive of phenomena that arise when a symmetry is broken [371. 

In the example discussed above, we took for the internal Hamiltonian Eq. ( 
Hamiltonian may be engineered by "dressing" a 2-level atom with a laser [T^llll]. Alternatively, one can dress a spin- 
1/2 system that has a magnetic moment with an external magnetic field, as is the case with Berry's Hamiltonian|16j . 
Below we describe a 2D version of such systems. 



where U is given by Eq. ( 27 ) . Such a 



E. Aharonov-Bohm gauge potentials 

The systems in the previous sections provides illustrations of how effective gauge forces manifest in simple quantum 
systems, and serves as a template for more intricate systems. Below we describe a somewhat more complex 2D system 



that has the identical structure to that described above, except that U given by Eq. ( 27 ) is now replaced by 



(p 11 <p 

U = exp(-2-cr3) exp(2-cri) exp(z-CT3) (96) 



where cj) the azimuthal angle in the x y plane. Hamiltonian Eq. (|4| is now given by 

H = 1- Had 

Zm 

Had = Ui^)HBo{p)U\^). (97) 

This adiabatic Hamiltonian is identical to that of a spin 1/2 particle, with a magnetic moment /io, interacting with 
an external magnetic field in xy plane as shown in Figure ([5|, i.e. 

Had = fJ-o^- ■ B 

B = 4>Bip). (98) 
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The strength of field -B(p), where p is the distance from the origin, is set to have the constant value A, the BO energies 
are simply given by ±A, where we have set /Lto = 1, as in the example discussed in the sections above. The geometry 
of a physical realization for such a system is illustrated in Figure [5] In the adiabatic picture we again obtain Eq. (21 ) 
but with the vector potential now given by 



2p 



I e 



-i4> 



1 



(99) 



where (j) is the polar angle in a cylinder coordinate system. The diagonal scalar potential V = Hbo contains the 
BO eigenvalues as its entries. Eq. (99) contains diagonal as well as the off-diagonal components that couple the two 
adiabatic channels. Suppose the channel with BO eigenvalue +B{p) = A is closed, i.e. the collision energy fc^ < 
as p -> oo. 




FIG. 5. (Color online) Illustration for the system described by Eq. (971. Dotted line represents a classical trajectory of the 
neutral atom for a given impact parameter. Green concentric circles represent the external magnetic field lines. Blue arrows 
represent the spin orientation of the lowest energy Zeeman level for the atom. The upper Zeeman level is closed, and so the 
spin of the atom is slaved to its lowest Zeeman level. 

We then ignore the gauge coupling to the closed channel and project onto the open sector to obtain the approximate, 
one-channel, equation 



2m 



(V - iA)^F{R) + V{R)F{R) = EF{R), 



(100) 



17 



where 

2p 

ViR) = -B(p) + ^. (101) 

This procedure is analogous to the one we carried out in the previous section for the 2D "magnetic" slab. There we 
showed that the BO projection is accurate in cases where the excited channel is closed, this conclusion also applies 



to this system ^38j. Vector potential (101 1 has the property 

dr-A^O (102) 



c 



where the line integral is taken along the path C that encloses the origin. Because of Eq. ( 102 1 we cannot find a 
gauge transformation in which the vector potential is eliminated [26]. Indeed this vector potential is familiar from 
Aharonov-Bohm scattering [321 . A similar vector potential has been invoked for poly-atomic systems in which 
conical intersections occur. It is note-worthy to mention that the adiabatic eigenstates U{9)\± > (where |± > are 
the eigenstates of Hbo) are not multi-valued, and the BO eigenvalues zLB{p) are non-degenerate everywhere, so the 
analysis given in Ref. |151 140j does not apply here. Previous studies of geometric phase effects in reactive scattering 
require the presence of conical intersections, as in studies of the hydrogen-exchange reaction pT]. 



Though vector potential Eq. ( 101 1 is nontrivial it does not lead to an effective Lorentz force for p ^ 0, nevertheless 
its presence has a profound effect on the scattering properties of this system[38,. These examples illustrate how gauge 
potentials arise due to the "dressing" of the atom (a spin 1/2 system) by external fields. Below we show how this 
relates to a central question of geometric magnetism, how and why do non-trivial gauge potentials arise in inter- 
atomic interactions? We now take our system to be that of two spin - 1/2 systems that interact via a spin-spin dipolar 
interaction. Such systems often realized in nature, including that for the interaction of two hydrogen atoms [201 138) . 
They also serve as a template for a non-trivial two-qubit quantum gate [20]. Here the adiabatic energy splittings are 
generated by the internal spin-spin interaction, and the BO energy matrix is sandwiched by rotation operators similar 
to that given above. In an adiabatic expansion (now a four channel, or two-qubit problem) the PSS equations possess 
a four channel vector potential, similar to those given above and discussed in detail below. 



IV. INTERNAL GAUGE POTENTIALS 



Consider Hamiltonian ( [98^ in three spatial dimensions and with B spherically symmetric about the origin. This 
(or closely related to it) system has been discussed extensively in the literature Jj6, ,23, j42j ■ In Berry's model ^16^, R 
is treated as a classical variable, Stone quantized R and considered the full quantal solutions. Berry and Robbins]^ 
also studied, from a semi-classical perspective, the role of the effective Lorentz forces[ll[51 that arise in such systems. 
The diagonal components of A produce the effective magnetic field of a Dirac magnetic monopole[2, 43^. These gauge 
potentials are realized in molecular and atomic systems [21 [18] and Stone ^Jj suggested a "Heath- Robinson- type 
device", a rotating solenoid interacting with a spin 1/2 particle, to serve as a physical realization. 

Consider the Hamiltonian for two, interacting, spin- 1/2 particles [20] with magnetic moments 

H=^+ HadiR) (103) 

where R = {R, 9, cj)) is the 3D inter-atom separation vector expressed in polar coordinates, and p is the relative 
momentum for the two particles. The adiabatic Hamiltonian is given by |44j 



Had — Ho + Hdip 

Ho^ (^j:{R)~^j:{R))Sa- 



•>b 



33S(i?) +1 E(i?) 



Hdip - ^3 



Sa ■ Sb — 3 



4 

{Sa ■ R){Sb ■ R)- 



R^ 



(104) 



The scalar potentials, ^^{R), ^S(i?), are arbitrary and depend on the polarization properties of the constituents. 
For a structure-less system, e.g. two electrons, these scalar terms are degenerate and Hq reduces to the long-range 
Coulomb interaction. In the case, of two neutral spin 1/2 atoms, such as two Alkali atoms, Hq (in atomic units) 
represents the isotropic electrostatic interaction between the atoms and ^Y,{R) and ^S(-R) are the triplet and singlet 
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BO energies of the electrostatic Hamiltonian respectively. Si are spin-1/2 operators. Hdip is the long range component 
of the anisotropic magnetic interaction, a is a constant and for two electrons is given by the fine structure constant 
when atomic units are employed. We re-express [201 Hi] 

Had = U{e^)HBoiR)Uie(b)^ (105) 

where 



U = 

Ua 



and 



H 



BO 



\ 



2i?,3 



'm) - 






/ 



(107) 



is a function of the internuclear distance R only. Using definition ( 11 ) we find that the induced vector potential has 
the block form, 



A = 



At 




(108) 



R 



tan^^ 



CXp(— 10) 

V2 



exp(i0) / 
2J 72 

{ie + 0) 







CXp( — 

%/2 



ie + 4)) \ 



■{ie + ct)) -</)tan[|] 



We note the following features [201 H5] 



C/fe = cxp(— icr30) exp(— i(T2^) exp(itT30) 
/O 1 \ 

-75 -75 
-71 71 

V 1 / 



(106) 



• The degenerate triplet sector of the diagonal scalar potential is split due to the inclusion of Hdip. The states 
corresponding to magnetic quantum numbers \M\ — 1 contain an attractive component —cP'IIR?^ whereas the 
M = state contains a repulsive component o? j^R? . 

• Non-trivial vector gauge potentials are induced. They contain off-diagonal components that couple the M = ±1 
adiabatic channels with the M = channel. There is no direct coupling between the two M — ±1 channels, 
and the diagonal components are that of Dirac-magnetic monopoles of integer charges ±1. 

• Because the energies of the diagonal entries in Hbo are split, a gauge transformation cannot be found in which 
the temporal components of the field strength tensor vanish 0]. As a consequence the gauge potentials are 
non-trivial in the sense described for the systems discussed in the previous sections. 

When solving the scattering equations in the adiabatic gauge the induced gauge potentials, depending on the collision 
energy, contribute to spin changing cross sections. Though exact scattering calculations for such a system (i.e. two 
ground state hydrogen atoms) have been performed in the diabatic gauge ^45) the question of how and in what manner 
the induced gauge potentials infiuence scattering cross sections remains to be investigated [20 j. 



V. SUMMARY AND DISCUSSION 



In molecular physics and atomic collision theory one is typically faced with a many-body Hamiltonian having the 
structure given by Eq. Q. The adiabatic Hamiltonian iJa(i(-R, f") is a function of the "fast degrees" of freedom r 
as well as i2, the coordinates for the positions of the atoms. Using the eigenstates of Had as a basis one arrives at 
the set of Equations (|3]) in which the matrix vector gauge potential A manifests. In applications, Eq. ([s]) involves 
an infinite set of coupled equations and a truncation procedure must be invoked in order to proceed with a working 
theory. The remaining gauge potentials, after truncation, are non-trivial in the sense that the resulting spatial field 
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strength tensor does not vanish[l] and, as a consequence, lead to effective gauge forces. They emerge in addition to 
the traditional Hellmann-Feynman force. We introduced and studied the gauge structure of several low-dimensional 
scattering systems that are solvable and where it was not necessary to employ the truncation, or BO, approximation. 
For the adiabatic Hamiltonian we replaced the fast degrees of freedom with the discrete Pauli coordinates cr, 

HadiR,r) ^ HadiR,(T) 

so that a complete set of basis functions can be used to obtain the exact solutions of the PSS equations (|3|. Though 
our choices for the adiabatic Hamiltonians varied in complexity they shared the feature that 

Had = U{R)HBoU\R) (109) 

where U is a unitary operator which can be characterized by it's gauge field as given in Eq. ([6|. Hbo is a diagonal 
matrix that contains the energy eigenvalues for Had, and in our discussion we were careful to choose systems in which 
Hbo did not feature degeneracies at physical values for the slow parameter R. The energy gap in Hbo explictly 
"breaks" the symmetry inherent in the definition of U . The matrix structure of Had depends on the particular system. 
For example, in our discussion relating to Figure ([2| Had can be "engineered" [T^] by the application of a laser field 
on a two-level system, for the system described in Figure ([5| Had is generated by an external magnetic field coupled 



(104 1 for two spin-1/2 particles. 



to a neutral spin-1/2 particle. In Eq. (106) [/ is a rotation operator that characterizes the dipolar interaction Eq. 



For the system described by Eq. ( 28 1 we were able to solve for the scattering parameters in both the adiabatic and 



diabatic descriptions and showed that the quantity, 

F= I ds-j (110) 



c 



where j is a current, is invariant with respect to the choice in representation. Furthermore, for the closed loop C shown 
in Figure (|4| we found that J- does not vanish in the low-energy limit. In the adiabatic representation the non- vanishing 
contribution arises from the gauge potential A that define U. In the low-energy limit, where the excited adiabatic 
channel is closed, we found that the BO projection procedure serves as a good approximation for solutions obtained 
by the fully coupled calculation. The efltective single-channel gauge potential predicts a "magnetic" field, given by Eq. 



(38), that leads to the appearance of an effective Lorentz force on the scattered atom. The BO amplitude which is 
coupled to the, non-trivial, gauge potential allows the BO current to satisfy the condition T ^ Q. We compared the 
quantal scattering solutions of this system with that of a charged particle interacting with a ferromagnetic slab, and 



we found that the former predicts, according to Eq. (46), a deflection angle in harmony with that predicted by the 
latter. Though the multichannel gauge potential A is trivial in the sense that the spatial components of the curvature 
vanish, in accordance with Eq. ([8|, non-trivial gauge behavior manifests in the low energy limit. The adiabatic 
Hamiltonian does not suffer from any degeneracies, and we propose that the non-trivial gauge behavior is akin to 
the phenomenon of symmetry breaking 13 7|. This system serves as a valuable template for more complex examples, 
but it also highlights possible applications for geometric magnetism. For example, because of the equivalence for the 



predicted deflection angle Eq. ( 46 ) with that of scattering by ferromagnetic strip we can anticipate its application as 
an effective "magnetic" lens for neutral atoms [46]. 

We showed the relationship between systems in which dressing is carried out by external agents, such as a magentic 



field for the system described by Eq. (97), or a laser field Jj] for the two level system given by Eq. (27), with those 



in which it is accomplished by mutual inter-particle interactions. For the interacting two-qubit system give by Eqs. 



(105) and (106) "dressing" is accomplished by the 2-body tensor interaction Eq. (104). We argue that, regardless of 
the nature of the dressing agent, gauge behavior manifest in these disparate systems via the same mechanism, and 
we do not advocate the use different labels, e.g. fictitious, artificial, synthetic or geometric magnetism, to identify the 
phenomena. We propose that geometric magnetism is an appropriate terminology that underscores its universality. 



Appendix A: Coupled channel solution 



We compare the solution of the approximate adiabatic solution Eq 
equations (21). For the sake of economy in notation we set h = 1 in this appendix 
ansatz 



( p4| to that obtained by solving the coupled 
In region a; > we take the 



exp(zri x) 



(Al) 
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Inserting this ansatz into Eq. (22 1 we obtain a set of homogeneous hnear equations for the coefficients CotCc whose 



non-trivial solutions require that Q. obey the eigenvalue equations 

-VL^ - Al- Al + 2mE =F ^Ji2mA - 2Aon)^ - Ajfl^ = 0. 

Analytic expressions for the eigenvalues can be gleaned using the formula of FerrariajlJ for the roots of the quartic 
equation. A much simpler and transparent expression that is valid in the limit where m — > oo is sufficient for our 
purposes. In this limit the roots are given by 



An±i\ 2m\A - E\ 



Ai 



(1 + S/A) 



-^0 ± A/2m(i? + A) - ^(1 - E/A). 



(A2) 



2rn 



and require that << A i.e. the collision energy is much smaller than the energy defect 



We set E 

between the two BO potentials. We then obtain the following solutions to Eq. (|21|, for a; > 0, 

1 / An 



= eyiY>{iA(3x) exp(±V 4mA : 



.JAI + AI+ 4Am 



- 2i\/ Am 

A2 



and 



F_g = exp(— i^lo^;) exp(±i J — A\ 



In the region x < Q we obtain 



F = 



Sco exp(\/4mA x) 
cxp(ifc(a; — L)) + Re:sq){—ik{x — L)) 



(A3) 



(A4) 



(A5) 



We solve for the open channel reflection coefficient R by requiring the solution at a; > to match the boundary 
condition at x — L, and then matching it with the asymptotic solutions at a: = 0. However, care must be taken 
when matching the derivatives of the adiabatic amplitudes since the vector potential A is discontinuous across the 
boundary. We require the diabatic amplitudes given in Eq. ( 18 1 and their derivatives to be continuous across the 



boundary. However since the adiabatic amplitudes are related to them according to Eq. ( 20 1 we require that 



F' 



F' 



(A6) 



where F'^ady^ £^'ad< derivatives of the amplitudes in the regions a: = ± e respectively, and F^adi^) ^^'^ value 

of the amplitude at a; = 0. Matching the amplitudes and enforcing condition Eq. (A6) we obtain 



R=-l + 2ik{L- 



tanh(AiL) 
A, 



(A7) 



Appendix B: Scattering by a 2D ferromagnetic strip 

We consider the scattering of a charged particle (here we set ^ = 1) by a 2D ferromagnetic slab as shown in Figure 
[6j The particle with energy E is traveling to the right, along the a;-direction, and impacts the slab, in the normal 
direction, of thickness L. We consider the Schrodinger equation 



2 m 



The vector potential describing system is given by 

A = j \Ba9{x + L/2){x + L/2)e{L/2 - x) + e{x - L/2)LBa 
■id Ay 



B 



OA, 



zBoe{x + L/2)e{L/2~x). 



\ dx dy 

We assume the particle approaches from the left (for x < —L/2 ) with the incident wave boundary condition, 

Tp = exp(ifca;) + r exjp{—ikx) 



(Bl) 



(B2) 



(B3) 
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where k = ^j^E and r is the reflection coefficient. In the region — L/2 < x < L/2 the Schrodinger equation 
becomes, 



whose solutions are, 

ciD^^^ (y/Bo/2{L + 2x))+ csL* .2+^ (iy/Bo/2{L + 2x)). 

2Bq ^ / 2Bq V / 

where Djj{x) are the Parabohc Cyhnder Functions ^5]. For x > L/2 the ampHtudes satisfy, 

_ + fc2^ = 

ox'^ 

and whose form is determined by the scattering boundary conditions 



^(x) = t exp(i V fc^ - ^'^x) fc > |$| 
i^ix) = t exp(-v^|fc2 - $2|2,) fc < |$| 



1. Reflection and transmission coefficients 



(B4) 
(B5) 



(B6) 



(B7) 



We now solve for reflection and transmission coefficients. We take the logarithmic derivative, at a; = L/2, of the 
internal solution Eq. (B5| and match it with that obtained using the asymptotic boundary condition Eq. (B6) 
resulting in the ratio. 



D'. 



/2B^L) -iy/W^^D 



rlB^L 



7 = C2/C1 = 



-D'_J^_^ {^/2B~qL)+i./W^WD_j^_, {V^L)' 

2Bq 2 2Bq 2 

Taking the incident wave and matching it, and its derivative, with the internal solution we obtain 

k — iy 

2Bn 2 2Br, 2 

y 



r = exp{—iLk)-_ 



D 



7D 



and 



t = 2jA:exp(— j — ) exp(— — V k^ — $2) x 

(d^2_s, {V2BoL)+jD ,2^s„ {iVm)L] 

(Dj^ 1 +jD fc2 

\ 2Bq 2 2Bg 



D\2 



ik 



where we used the shorthand = D^{Q). 

Consider the limit L — > as the flux density $ = BqL stays constant. We flnd, in this limit, 

as 



y — ^ fc2 - 



Thus 



k-y/k^- $2 

k + Vk'^- $2 

t = l + r 



(B8) 



(B9) 



(BIO) 



(Bll) 
(B12) 

(B13) 



and are the scattering parameters obtained for a system in which a particle is incident on a step function potential of 
height $2 at X > 0. 
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2. Conserved currents 



For the Schrodinger equation 



the current density is given by 



'^9^ = -2^r-^^) 



ih 

j = (iPVij* - ^*Vip + 2iAi}*i;). 

2m 



Inserting expression Eq. (B2) into Eq. (B15I we obtain 

hk 



Jx 



m 



which can be expressed as an incident flux ji 



and a reflected flux jr 

h 



m 



3v 



(B14) 
(B15) 
(B16) 

'|^|rp. In the region x > L/2 we find 

(B17) 



where we have used Eqs. (B2|, (B7) and (B16). Thus the transmitted flux makes an angle 9 with respect to the 
normal given 



tan^ : 



Jy 



where — BqL, and agrees with the classical angle of deflection. 



(B18) 
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FIG. 6. (Color online) Normal incidence of a charged particle on a ferromagnetic slab of width L. The magnetic field B is 
directed into the page and the dashed line depicts a typical classical trajectory. In this figure the classical trajectory describes 
a particle that does not have the required velocity to overcome the barrier shown in Fig. |2] 
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